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ABSTRACT 



It is desirable to perform a study of the elastic 
responses of submarine hulls and internally attached 
structures to underwater shock loading. In view of the 
expense and possible involvement of non-shock testable 
equipment it is not always desirable to perform actual shock 
tests utilizing real hulls or SSTV's (Submerged Shock Test 
Vehicles). With the advent of large scale computing power, 
numerical methods now exist to predict equipmen t/hul 1 
responses. The ELSHOK, (ELASTIC SHOC K) computer code is 
used to perform a parametric study on a submerged shell with 
an internally att ached . substructure . Of particular interest 
is the phenomenon of dynamic amplification or lack of such 
in the response of this particular substructure and effects 
on substructure-shell interaction due to varying 
substructure mass and stiffness subjected to underwater 
explosion side-on shock loading. 
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INTRODUCTION 



A. BACKGROUND 

A submarine must be able to withstand moderate to severe 
shock loadings that will result from und erwater explosions 
(UNDEX) it might be subjected to in modern warfare. 
Conventional and nuclear weaponry deliver devastating forces 
and no vessel can be expected to survive in the events of 
near or direct hits. The emphasis on design and testing is 
in the area of survival of mission critical equipment, 
machinery and weapons systems of a platform. The ability of 
a naval vessel to carry out its mission after being 
subjected to an UNDEX will depend on the survivability of 
these installations. 

The current specifications for building ships and 
submarines contain the requirements for shock loading to be 
met by the builder or vendor of shipboard installations 
and equipment. In general, all critical equipment is 
required to pass a series of shock tests where testing is 
practical as outlined in MIL-S-901D fPef. 11. This document 
specifies the shock testing requirements for shipboard 
machinery, equipment and systems which must resist High 
Impact (HI) mechanical shock. The purpose of these tests is 
to determine the suitability of machinery, equipment and 
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systems for use after exposure to severe shock which may be 



incurred in wartime. It is not always practical -to test 
equipment due to size and weight. In the case of equipment 
design, it is not always practical to evaluate an 
installation prior to final design. Mechanical shock 
testing of individual equipments may not always indicate how 
equipment will perform or respond when actually installed. 

An analysis for a given design is specified for the 
situations where testing cannot be conducted. It is based 
on limited and dated testing and observations which are used 
to specify a shock desigin factor. A static design force is 
arrived at for various weights of equipment which is used to 
evaluate the design of support structures of installations. 
The method is simple to apply but does not always correlate 
well with actual testing. 

Numerical methods have been developed for use in 
submarine shock response analysis. The computer code ELSHOK 
(Elastic SHo ck ) developed by Weidlinger Associates for the 
Office of Naval Research [Ref. 2] is used to conduct this 
study . 



B. PURPOSE 

ELSHOK has been developed through a series of controlled 

tests and has been used to predict equipment responses prior 
to low level explosive tests on submarines. Although ELSHOK 

has been validated by extensive testing it is considered 
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prudent to conduct additional research to gain insight to 
the sensitivity of the code's predictions when a variety of 
conditions are applied. The intent of this investigation is 
to examine the .change in the equipment/hull interactions 
under these conditions. The models used to conduct the 
parametric study are similar to those used in the mid 1970's 
in tests conducted by the Office of Naval Research (ONR) 

[Ref. 3]. The internal structure used was specifically 
designed to reduce the effects of dynamic amplification for 
the conditions applied in these tests. The same model 
geometry was studied here by varying one main dimension to 
tune the substructure to the natural frequencies of hull to 
see if this particular internal structure configuration 
retained this design feature over a range of mass and 
stiffness variations. 
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UNDEX SHOCK RESPONSE USING ELSHOK 



1 1 . 

A. ELSHOK PRINCIPLES OF OPERATION 

The ELSHOK computer code explained in [Ref. 2] is a 
suite of computer programs assembled and developed to 
calculate the submerged shock response of a linearly 
elastic structure subjected to an underwater shock wave. 
Specifically, a ring-stiffened shell of revolution of finite 
length, with or without internal structure, is considered to 
be immersed, initially at rest, in an infinite acoustic 
fluid and to be excited by an acoustic shock wave. The 
shock wave originates from an UNDEX arbitrarily located from 
the structure and is assumed to propagate as a spherical 

r 

wave through the medium. This assumption allows inciden.t 
fluid particle velocities impinging on the wet surface of 
the immersed structure, the hull of a submarine, to be 
determined from the incident fluid pressures. ELSHOK 
contains two possible representations of incident pressure. 
One uses a known pressur e-t ime profile. This type is of the 
kind that may be used to describe a tapered charge simu- 
lating a nuclear explosion. The other is the empirical 
decaying exponential put forth by Cole [Ref. 4] for spheri- 
cal charges of conventional explosives. If the shell 
contains internal equipment, a substructuring technique 
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described by Ranlet [Ref. 5] treats the internal equipment 
response by coupling the free-free modes of the empty ring- 
stiffened shell and the fixed-base modes of the internal 
equipment through use of dynamic boundary conditions. 
Structure-fluid interaction is calculated using the Doubly 
Asymptotic Approximation (DAA) presented by Geers [Ref. 6]. 
Normal fluid displacement of the structure-f luid interface 
is expanded in a series of surface expansion functions which 
are orthogonal over the entire wet surface of the submerged 
hull or shell model. The surface expansion functions lead 
to matrices in which the elements are determined by matching 
the limits of zero and infinite frequencies of the pressure- 
velocity relationships. Thus, exact solutions are obtained 
for the transient problem at early and late times and DAA 
results in a smooth transition between the two limits. 
Because DAA calculates the effects of the fluid in the 
interaction problem through the series of surface expansion 
functions defined only at the wet surface, the fluid field 
is in effect uncoupled from the structural field. The fluid 
effects produced by DAA are simply applied as additional 
loading in the modal component structural analysis of the 
shell . 

The substructuring technique separately determines the 
vibration modes of the components that comprise the combined 
submerged structure of shell and internal substructures. 
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Equations of motion of the entire structure are obtained by 
combining interaction forces and moments and enforcing 
compatibility of deformation at component connection points. 
This precludes calculations involving the modes and natural 
frequencies of a system of combined components and also 
avoids a requirement for a combined system stiffness matrix. 
The ELSHOK code also handles the case of no substructure in 
the shell using the same substructuring method. 

ELSHOK is particularly suited to the task of a 
parametric study due to its purely numerical nature and 
modular organization. Various internal equipment models can 
be combined repeatedly with the same shell model. Inputs of 
pressure- time profiles or decaying exponential types are 
mathematically applied to these numerical models and result 
in output in the form of veloc i ty-t ime histories. This 
transient response analysis takes into account interactions 
of shell and substructure and eliminates dependency on data 
derived from actual explosion testing. 

B. ORGANIZATION AND IMPLEMENTATION OF ELSHOK 

The implementation of ELSHOK is accomplished using a 
suite of seven major programs. These major components are: 

1) B0S0R4 - structural analyzer for shell [Ref. 7] 

2) ACESNID - virtual mass processor 

3) PIFLASH - shell-fluid processor 

4) SAP IV - structural analyzer for substructure [Ref. 8] 
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5) PICRUST - substructure processor 

6) USLOB - time integration processor 

7) PUSLOB - plotting processor 

Illustrated in Figure 1 are the general relationships 
between the main components and four phases into which 
ELSHOK is separated. In any analysis using ELSHOK the shell 
model and internal structure model must be available or 
constructed by the user. Development of the models used in 
this study is discussed later. Once models are arrived at 
the analysis can be carried out in the four phases by 
executing each of the seven major programs sequentially in 
the order listed above. 

1 . Phase I--Shell and Fluid Analysis 

The finite element difference code B0S.0R4 is used to 
analyze the idealized model of the submerged hull or shell. 
The model must be symmetric about its longitudinal axis and 
it may contain rings and bulkheads. Circumferential rings 
can be treated as discrete stiffeners or, if closely spaced, 
they may be represented in the orthotropic approximation 
(smeared) [Ref. 9]. An entire main body or full model is 
used to capture the gross effects of the shell response when 
calculations are made for circumferential wave numbers N = 0 
and 1. These include rigid body translation and whipping of 
the structure (N = 1) and the torsional modes of the struc- 
ture (N = 0). Performing the analysis for N = 0 and > 2 
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(breathing mode N = 0) obtains the local shell response. In 
this analysis wave numbers N = 0, 1, 2, and 3 are used. 
Figure 2 shows the basic nodal patterns for N = 1-4. ELSHOK 
requires separate runs of BOSOR4 for each wave N used in the 
analysis. Each BOSOR4 run determines the in-vacuo free-free 
modes and corresponding natural frequencies from the dynamic 
properties and geometry supplied in the model of the shell. 
If the shell is symmetric about the midpoint of its length 
and the shock wave impinges on it so as to produce symmetric 
reponse, only half of the structure need be modeled. Such 
is the case in this study since the shock wave emanates from 
an origin located normal to the shell at midships. 

ACESNID (Accession to I_nertia and Damping) accom- 
plishes the second step of Phase I by calculating the 
virtual mass array. This provides the late-time contri- 
bution of the Doubly Assymptotic Approximation [Ref. 6]. 
Normal displacements corresponding to the surface expansion 
functions are applied to the surface of a cavity in the 
fluid (infinite medium) having the same shape and size as 
the wet surface of the shell model being investigated. The 
cavity geometry is obtained from the shell geometry 
contained in the BOSOR4 output files. Only one run of 
ACESNID need be performed to encompass all values of N being 
used and all values of M being used in BOSOR4 must be 
indicated in ACESNID for compatiblity . 
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PIFLASH (Prepare Input for Fluid _and Shell ) 
accomplishes the last step of Phase I and that is to 
reorganize the fluid and shell data files from B0S0R4 and 
ACESNID into a fluid-shell file that is used to solve the 
response equations. If there are no internal substructures 
involved in the analysis, then this fluid-shell file is the 
main input for Phase III and Phase II is not necessary. 

2 . Phase II — Structural Analysis of Internal 

Subs true tures 

A finite element model for each internal 
substructure must be developed. The SAP IV code performs 
the structural modal analysis on each substructure to be 
included. The code contains a variety of types of elements 
that can be used with ELSHOK. Concentrated masses attached 
to the substructure may also be represented by the methods 
used in SAP IV. Each substructure requires a separate 
performance of calculations in which the fixed-base mode 
shapes and natural frequencies for each are determined. 

This data is recorded in the "substructure mode file" along 
with the geometry, stiffness and mass information of the 
model for future processing. 

PICRUST (Prepare I^nput for Calculating Response of 
Substructure) rearranges the data from the "substructure 
mode file" and computes additional information required to 
solve the response equations. PICRUST also provides for the 
selection of which modes are to be included from SAP IV and 
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the specifics of connectivity between the shell and 
substructure used in the response calculations. 

Additionally, influence coefficients corresponding to the 
base motions of the substructure, constraint inodes and modal 
coefficients for calculation of interaction forces with 
moments between the shell and substructure are determined 
[Ref. 9]. PICRUST is the last step in Phase III and its 
results are stored in the substructure file. 

3 . Phase III — Shock Response and Time Integration 
The information required to complete the fluid- 
structure interaction is supplied as one input file. 
Therefore, if a substructure is being analyzed in the 
problem the shell-fluid file from Phase I must be combined 
with the substructure files produced in Phase II and become 
the input file for USLOB (Underwater Shock Loading of 
Bodies) which is ELSHOK's time integration processor. For 
the case of an empty shell the shell-fluid file is the input 
for USLOB. Integration is accomplished using a modified 
Runge-Kutta scheme. The time step and shock wave loading is 
specified by the user in USLOB. Care must be taken to be 
sure that the product of the time step chosen and the 
highest circular structural frequency w of the modes being 
used from SAP IV is less than unity. Otherwise, the 
numerical integration scheme will not converge without 
changing the time step. 
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The cases studied in this work used both methods 
available to represent the incident fluid pressure. A 
pressure-time history profile shown in Figure 3 was supplied 
as a typical tapered charge simulation. Nine discrete 
pressure- time points are entered as input to describe the 
incident shock wave loading. Tapered charges are used in 
UNDEX testing to produce the same type of loading resulting 
from an underwter nuclear detonation. Here the impulse of 
the shock wave is delivered to the target over a longer 
period of time compared to the interval over which the 
impact of a conventional explosive charge is delivered. 

Both are assumed to emanate from an origin far enough away 
so as not to include the target in the gas bubble effects. 
ELSHOK is used for elastic response analysis only and hulls 
within the immediate vicinity of the gas bubble formed in an 
UNDEX would most likely have responses well into the plastic 
response range. Figure 4 shows the pressure- time profile 
used to represent the incident shock loading of a spherical 
charge of HBX-1, a conventional explosive, of a weight which 
would deliver the same impulse over the same period of time 
as the duration of the tapered charge loading and detonating 
at the same range from the closest point of the target or 
standoff. This profile is an empirical exponentially 
decaying pr essur e-t ime history where the initial pressure is 
the peak incident pressure on the shell due to the charge 
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weight and standoff. It is seen from Figure 4 that the 
majority of impulse of the HBX-1 blast is delivered in'side 
the first millisecond. The empirical equation over which 
USLOB performs integration to obtain the fluid particle 
velocities is of the form [Ref. 4]: 



Pj(R,t) = K 1 (W X/J /R) exp ( - t/ 6 q ) 



( 1 ) 



where, 

P j (R, t) = the incident pressure on the hull at range R (ft) 
from explosive source to point of interest (psi) 

= multiplicative constant for incident pressure 

= spatial decay constant for incident pressure 

= time after arrival of shock wave at point of 
interest (msec) 

= weight of spherical charge (lb.) 

K 



K- 



K. 



W 



K. 



K 



= (W^^^/R) ^ , time constant of 

exponential decay (msec) 

= multiplicative constant for time constant of 
exponential decay 

= spatial decay constant for time constant of 
exponential decay 



( 2 ) 



USLOB creates a generalized velocity file for the 
shell-substructure and shell-substructure transformation 
file to be used in the last phase of ELSHOK in producing 
plotted and printed velocity-time profiles. The trans- 
formation file contains information required to describe the 
motion at the shell-substructure connection points. 
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4. 



Phase IV — Plotting Velocity Profiles 



PUSLOB (Plots for USLOB ) processes the results of 
the previous three phases so they may be plotted or printed 
as velocity-time histories. Velocity time-profiles are 
requested for various locations on the shell and 
substructure by specifying a call out point in the input 
files. The call out points are the nodes that are used to 
describe the models of the shell and substructure in B0S0R4 
and SAP IV. Selection of nodes is important because 
velocities can only be obtained at those locations. The 
shell velocities may be represented in either a global or 
local coord inate . system. The global is represented by X, Y, 
Z that locate the positive sense of X toward the bow or 
forward section of the shell, the positive sense of Y is to 
port and positive Z is up. The local system is described by 
u-longitudinal (fwd pos) v-circumferential ; w-inward normal 
from the shell. The substructure is always described by the 
coordinate system specified in its SAP IV finite element 
model . 



21 



PHASE n 
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Figure 1. Organization of the ELSHOK Code [Ref. 2] 






Figure 2. Modal Patterns for Cylindrical Shells 
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Figure 4 



Pressure-Time History for the Conventional Charg 
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MODELS INVOLVED IN THE ANALYSIS 



A. MODELING WITH FINITE ELEMENT/FINITE DIFFERENCE METHODS 

Finite element/finite difference methods have become 

powerful tools in engineering with the advent of large scale 
computing power. These methods are of great use in modeling 
techniques. Even though the models are much simplified for 
use in these techniques the results arrived at can be quite 
accurate or at least adequate for the problems that might be 
impossible to solve by other methods. Accuracy of the 
results is closely tied to how well the model represents the 
actual configuration, so care must be taken when 
constructing the model. In ELSHOK, the finite difference 
code B0S0R4 is used to describe the dynamic behavior of 
underwater hulls by representing them as much simplified 
models of ring-stiffened shells of revolution. Similarly, 
the models used to analyze internal substructures are also 
simplified representations made with the finite element code 
SAP IV. 

B. SUBMERGED SHELL MODELING 

Figure 5 is a schematic of the general shell structure 
and substructure model depicting the notation and convention 
used to describe the shell global coordinate frame X,Y,Z, 
shell local coordinate frame u,v,w and the substructure 
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coordinate frame x,y,z. The z-axis orientation of the local 
substructure coordinate system is related to the Z-axis of 
the shell by the angle where a represents the substruc- 
ture. Points are located along the meridian of the shell 
using the variable s. The angle 0 is used to locate 
circumferentially the desired local shell displacements. 

One shell model is used in all calculations of this study. 
Once the symmetric full model was constructed no changes 
were required to conduct the investigation of varied 
substructure sizes and types of UNDEX. 

The shell or hull that has been modeled for analysis by 
ELSHOK in this study is similar to an actual structure 
configured for testing by the Office of Naval Research 
[Ref. 3]. Figure 6 is a schematic of the representation of 
the full shell model used. Though its size is much smaller 
it is also of similar geometry to that of SSTV's (Submerged 
Shock Test Vehicles) that are used to investigate the 
underwater shock effects on submarine installed equipment. 
The model is an idealized submarine hull in that both can be 
considered as free-free ring-stiffened cylinders with end 
closures. The shell in this study is in fact a cylinder of 
revolution made of high strength steel stiffened by six 
large steel rings plus forty small stringers also of steel. 
The end plates are aluminum. An internal substructure, also 
similar to one involved in the testing by ONR is mounted at 



26 



midships on the port side of the hull. In B0S0R4 the model 
is divided into segments of consistent properties contained 
in each cross section. When the structure is symmetric 
about the midpoint of its longitudinal axis only one half of 
the model need be configured. This is accomplished by 
specifying the appropriate boundary condition at the point 
of symmetry. In this problem symmetry is employed allowing 
the shell model to be divided into two segments. Segment 
one is the aft aluminum end plate and has nine nodes 
specified on its surface. Node 1/01, node one on segment 
one, is the center of the end plate. The other nodes move 
as concentric circles to the outer circumference. Segment 
two is the cylinder from the end plate to the midships 
point. Twenty-seven nodes are specified for segment two 
that are actually bands around its circumference. Node 
numbers increase from aft to forward such that node 2/27 is 
the midpoint on the cylinder. 

Concentrated masses and substructures internal to the 
shell must be modeled separately if they contribute signi- 
ficantly to the stiffening of the shell. Therefore, the 
large rings have been modeled as discrete stiffeners and the 
internal substructure is modeled separately using SAP IV. 

If there were no interest in the substructure response it 
could be modeled as a concentrated mass attached to the 
shell. There is only one substructure^ doubly overhung 
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cantilever beam, used in this study and there are no 
concentrated masses. The smaller stringers do not 
contribute as significantly to the response of the shell as 
do the large rings so they are smeared onto the shell to 
account for their mass [Ref. 9]. 

The cylindrical section of the shell is a surface of 
revolution having consistent properties throughout its 
entire length. This allows information regarding Young's 
modulus, Poisson's ratio, mass density, thermal expansion 
coefficients, plating thickness and whether or not rings are 
smeared to be specified at only one point for the segment in 
B0S0R4. Similar information is specified for segment one, 
the end plate. Tables I and II are listings of the natural 
frequencies and wave numbers for each mode calculated. 
Twenty-two N = 0 breathing modes, thirty N = 1 rigid body 
and whipping modes and thirty N = 2,3 modes were retained 
for analysis. The torsional modes for N = 0 were not 
retained to eliminate a rotational input at the base of the 
substructure model. Modal contributions corresponding to N 
= 2,3 normal, meridional and circumferential displacements 
were retained for this full model case. Figure 2 shows the 
basic shapes of the waves on the cylindrical cross section. 



C. INTERNAL SUBSTRUCTURE MODELING 

In ELSHOK calculations, a substructure model is formu- 
lated to determine the modes, masses and natural frequencies 
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of the internal equipment being analyzed. Phase II of 
ELSHOK utilizes the SAP IV finite element code to accomplish 
modeling of submarine-installed equipment. ELSHOK is 
written in program modules that execute independently and 
exchange data through the output files. This implementation 
lends itself well to the type of parametric study conducted 
in this thesis work. A variety of different internal models 
may be combined with the same shell model to gain insight on 
the interaction of the response between the shell and 
various substructure characteristics. Modeling the 
substructures in a finite element method more accurately 
predicts effects of the substructure' on the hull to which it 
is attached than if the model is represented as a 
concentrated mass. A concentrated mass model would not 
account for any substructure responses that might become 
dynamically amplified. 

The substructure modeled for this study is a doubly 
overhung steel beam mounted with steel brackets to the port 
side of the shell. Connection points are at the two large 
rings at frames 18 and 27, placing the beam symmetrically 
about the midships section. Only half of the beam structure 
requires modeling since symmetry may be employed. Symmetry 
is indicated in SAP IV by specifying the appropriate 
boundary condition at the point of symmetry on the model. 
This is similar to the procedure in B0S0R4 . The dimensions 
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and details of the model are shown in Figure 7. Figure 7 
also shows the orientation of the substructure within the 



shell and to the shockwave. The doubly overhung beam is 
represented using the three-dimensional beam element model 
of SAP IV depicted in Figure 8 which shows half the 
structure with the first mode for D = 4.2 in. superimposed 
on it. In SAP IV nodes are used to divide the model into 
separate elements and serve the purpose of call out points 
for veloci ty-t ime histories in Phases III and IV of ELSHOK. 
Modeling one half of the doubly overhung beam was 
accomplished using 21 nodes, 4 on the bracket and 17 on the - 
beam half length. The brackets are represented as purely 
extensional members (a spring with mass) and the intro- 
duction of a bending moment at the connection between the 
beam and bracket is eliminated. This is accomplished by 
specifying an end release node at node 12 that represents 
the connection point on the model. The interest in the 
response of the substructure is limited to the z-direction 
and restrictions on the degrees of freedom of the model were 
appropriately specified. This reduces the calculations 
required and produces responses in the direction of 
interest . 

In order to investigate the response interaction between 
the shell and substructures of various weights and charac- 
teristics the depth D of the beam is changed for .each case. 
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It should be noted that the width of the bracket varies with 
the depth of the beam by 1.5 x D. Developing the various 
sized models was accomplished by editing the SAP IV input 
code to reflect the new cross-sectional areas and moments of 
inertia for the beam and bracket elements. Other properties 
remained unchanged. 



D. SHELL/SUBSTRUCTURE RESPONSE MODELING 

The structural problem of ELSHOK is separated into two 
parts to avoid calculations involving the modes, natural 
frequencies and stiffness of a combined system of the shell 
and substructure. Even though separate modal analysis is 
carried out on each model in different phases of ELSHOK the 
responses determined for the shell and substructure are 
arrived at by coupling the free-free modes of the shell with 
the fixed-base natural modes of the substructure. In effect 
the shell and substructure can be approximately represented 
by a two-degrees of freedom mass-spring model as in Figure 
9. Here the mass of the shell is M2 and the stiffness 
characteristics are K2. The substructure model is 
represented by Ml and Kl. It has been shown for some 
situations using two single-degree of freedom models can 
produce nearly the same responses as a two-degrees of 
freedom model when the mass ratio is very small. This in 
general does not hold for systems of large mass ratios and 
is apt to produce marked differences in response between the 
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two methods if the same stiffness characteristics are used 
in both [Ref. 10]. ELSHOK results can be viewed with more 
confidence in predicting component responses because it 
utilizes the multi-degrees of freedom type model. It is 
preferred in this study presented due to the large mass 
ratios that are involved. 
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END CLOSURE 
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Figure 5. Typical Shell-Substructure Configuration [Ref. 2] 
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Figure 6. Ring-Stiffened Cylindrical Shell Model Used in Analysis 
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Figure 7. Substructure Model Used in Analysis 
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Figure 8. Substructure Model Depicted in SAP IV 
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Figure 9. Two-Degrees of Freedom System Model 



IV. ANALYSIS 



A. OBJECTIVES OF ANALYSIS 

The objectives of the analysis procedure are to 
investigate effects on component responses brought about by 
varying the mass and stiffness of the internal substructure. 
More specifically, it is divided into two major and one 
minor task. First of these is to determine how heavy the 
internal structure must be before it significantly affects 
the motion of the shell. Secondly, the dynamic amplifica- 
tion of the responses of the substructure and its effects on 
the shell are investigated. Lastly, it is of interest to 
know the magnitude of the forces generated in the support 
brackets. An integrating scheme is used to determine 
deflections and thereby estimate the forces. 

The shell and substructure are subjected to side-on 
shock loading from an UNDEX. The UNDEX is represented in 
two different manners. One is an incident pressure- time 
history simulating a tapered charge or nuclear detonation 
and the other is an empirical exponentially decaying 
incident pressure simulating a spherical conventional charge 
detonation. The standoff or range from the origin of the 
UNDEX to the closest point of the target is the same for 
both the tapered and conventional cases. The weight of the 
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conventional charge is selected to deliver the same impulse 
as that of the tapered charge. First, the transient 
responses of the empty shell are determined for a basis of 
comparison. Next, responses for a variety of internal 
structure sizes are calculated for both conventional and 
tapered charge cases. The results are produced in the form 
of veloc i ty-t ime histories and are obtained for points 
specified on the shell and substructure. 

B. ANALYSIS PROCEDURE 

The shell natural frequencies and modes for wave numbers 
N = 0, 1, 2, and 3 are found using the BOSOR4 code. The 

virtual mass array computed in ACESNID is done including 
each wave N specified in BOSOR4 for compatibility. In order 
to simplify the analysis and since the brackets are designed 
to resist movement out of the horizontal plane, the 
torsional and rolling modes are eliminated from considera- 
tion to prevent a rotational input to the base of the 
substructure. This is accomplished by setting the variable 
NTORSN = 0 and omitting any purely rolling or torsional 
modes for N = 0 from input in the PIFLASH code. All thirty 
modes for the other values of N are used. This shell model 
is used throughout the study in combination with each 
substructure. The BOSOR4 modes deleted and retained are 
indicated in Tables I and II. 

SAP IV was used to determine the fixed-base natural 
modes and frequencies as D was varied on the substructure 
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model. The investigation begins with a beam depth D of one 
inch and increases to twenty inches. It should be noted 
that twenty inches is nearly physically impossible to fit 
within the confines of the shell and was computed to 
investigate the effects of a large mass on the shell. Each 
SAP IV run produced ten fixed-base natural frequencies 
corresponding to different mode shapes of the doubly over- 
hung beam. Table III contains frequencies computed for the 
cases D = 1.0 f 4.2 and 9.0 in. Not all the modes determined 
in each case were used in the follow-on calculations because 
some of the high frequencies would not allow convergence of 
the integrating scheme used without changing the time step. 
All modes having frequencies of 3979 Hz and greater were not 
used in the analysis. These higher modes do not contribute 
significantly enough to warrant varying the time step to 
remain within the convergence limits. 

Results from each SAP IV run require processing by the 
PICRUST code for preparation to merge the shell and 
substructure data files with the fluid data files. The 
configuration of each substructure installation along with 
the modes to be used in Phase III calculations are 
specified. The substructure z-axis is related to the shell 
coordinate system by the angle shown in Figure 5. The 
meridional locations of the attachment points are also 
specified here. Each case investigated placed the 
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substructure attached to the port side or an angle of 180 
from the shell Z-axis orientation (exactly one half the 
distance between the top and bottom) . The brackets connect 
to the rings at frames 18 and 27 placing the substructure 
evenly about the midships section of the shell or node 2/27. 
Each case D is subjected to two types of UNDEX representa- 
tions having origins at 70 ft. (840 in.) from the shell 
abeam the starboard side at midships. The origin also is 
located in the shell's Z-X plane. After PICRUST prepares 
the shell/substructure-fluid data the time integration of 
the incident pressures experienced by the shell surface is 
performed to determine velocities. 

The USLOB code models the type of UNDEX by differently 
representing the incident pressures. The tapered charge 
pressure- time history is represented by inputting discrete 
pressure-time points that describe it. The conventional 
spherical charge is represented by equation (1) and 
specifying the appropriate constants for the type of 
explosive used. The weight of the charge specifies the size 
of the explosion and its resultant impulse and incident 
pressure on the shell. In determining a charge weight that 
would deliver the same impulse as the tapered charge 
pressure- time history equations from Cole [Ref. 4] and 
constants describing HBX-1 were used. Describing the 
impulse of unit area of the shock wave front up to a time t 
after its arrival at the shell wet surface by: 
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I (t) 



P (t) dt 



(3) 



$ 



where: P(t) = P Q exp(t/9 Q ) (4) 

Then integrating equation (2) over the time after arrival of 
the shock wave results in: 

1 (t) = Vo [1 " ex P(-t/V ^ (5) 

Substituting equation (2) for 0 Q and P Q where, 

1 /3 K 2 

P q = K 1 (W 1/J /R) (6) 

into equation (5) then setting it equal to the area under 
the curve in Figure 3 the charge weight arrived at through 
iteration is 352 lbm of HBX-1. The constants for HBX-1 
are determined from data recorded in feet and msec. This 
required and to be adjusted to represent the empirical 
relationships in inches and seconds in order to be con- 
sistent with the units used in this application of ELSHOK. 

The standoff and lateral placement of the origin of the 
UNDEX is specified in USLOB along with the duration of 
integration and the time step to be used. A time duration 
of 10 msec and a step of 3xl0 - ^ sec was used to capture the 
peak response in each case. Also no effects due to surface 
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cutoff or reflection of shock waves were considered. This 
situation represents a deep submerged test simulation where 
bulk cavition due to the surface doesn't occur and there is 
no bottom bounce or surface cutoff effects. 

C. PLOTTING PROCEDURE 

Plotting of velocity profiles is accomplished by PUSLOB 
in the last phase of ELSHOK. Each profile is a record of 
the velocity history for a single node specified on the 
shell or substructure. The plot is limited to two angular 
locations of the shell nodes and velocities are expressed in 
in/sec due to the units used in this study. The version of 
ELSHOK available at the Naval Postgraduate School is 
installed on the VAX/VMS system and plots were produced on 
Tektronix devices. In order to aid in comparing the 
velocities of different nodes the data files used to create 
the PUSLOB plots were processed using program PLOTCONV. 
Program PLOTCONV rearranges the data files and converts the 
units of velocity from in/sec. to ft/sec. This facilitates 
the use of plotting routines EASYPLOT or DISSPLA available 
on the school's IBM system and expresses the plots in the 
customary units of the shock community: velocity in ft/sec, 

time in msec. Program FISH1 is a program developed for use 
with the DISSPLA plotting routine. Copies of programs 
PLOTCONV and FISHl are found in Appendix B. 
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D. FORCES IN THE SUBSTRUCTURE SUPPORTS 

Since ELSHOK produces results in veloc i ty- t ime 
histories, force calculations require some post processing. 
Forces generated in the support brackets would be mainly due 
to extension and compression because the bracket design 
resists motion out of the horizontal plane. Due to this and 
to simplify the analysis, inertia of the brackets is 
neglected in the force calculations, though their mass is 
included in the shock response calculations of ELSHOK. 
Program FORCES found in Appendix C was used to convert the 
velocity histories to deflections using a Simpson's one- 
third integrating scheme. Uniaxial forces estimated are 
based on F = eEA where e is the estimated deflection of the 
bracket divided by the length between the connection points 
of bracket at the ring and beam junction. E is the modulus 
of elasticity of steel and A is the cross section area of 
the bracket, 1.5 x D. 
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MODES AND FREQUENCIES OF SHELL MODEL FOR WAVE NOS 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 



TABLE III 



SUBSTRUCTURE NATURAL FREQUENCIES 



Frequency (Hz) 



= 1.0 in . 


D = 4.2 in . 


n = 9.0 


25.46 


106.9 


228.7 


81 .63 


330.5 


616 .3 


337.4 


1354.0 


1728.0 


463.0 


1421 .0 


3032.0 


966.4 


2700.0 


4982.0 


1206.0 


4167.0 


6552 .0 


1802.0 


5475.0 


8951 .0 


2151 .0 


6728 .0 


9957.0 


2680.0 


8380.0 


12430.0 


3281 .0 


9939.0 


17810.0 
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V. RESULTS 



Numerous cases were explored in this parametric study in 
which the beam depth D varied over the range of one to 
twenty inches. ELSHOK generated a large amount of 
information about each case. The most significant results 
have been summarized in tabular and plotted forms which 
contain information on cases subjected to the tapered charge 
and the conventional charge UNDEX simulations. SAP IV 
predicted frequencies and modes for the tapered charge case 
with beam depth D = 4.0 in. that compared very closely with 
the ONR model previously tested. This provides a certain 
amount of confidence in the SAP IV model used throughout 
this investigation. 

A. SHELL AND SUBSTRUCTURE INTERACTION 

Even though the variation of substructure sizes included 
the cases resulting near the 600 Hz secondary frequency that 
was expected to induce dynamic amplification in the original 
ONR model there is no such response over the entire range of 
cases investigated. The case for D = 4.2 in. conventional 
charge resulted in a peak of velocities at the tip and 
midsection of the beam as shown in Figure 10 and Figure 11. 
However, this is not a resonant condition. This slight 
amplification can be traced back to shell modes 1, 2, 3 for 
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N = 2 and mode 1 for N = 3. Tables IV and V summarize the 
range of results for both charge types and, when viewed with 
Figures 12 and 13, there is some indication of transfer of 
energy between the substructures and the hull, but not much. 
There seems to be no significant interaction between the 
beam and shell. The basic shape of the velocity-time 
histories at the beam tip, midsection and junction remain 
nearly the same throughout the range of D studied. The only 
effects seem to be the decrease of beam and bracket maximum 
velocities as the weight of the beam is increased. The 
complete velocity profiles for the beam bracket and connec- 
tion points for the D = 4.2 in. cases are shown in Figures 
14 through 17. This size was selected for comparison 
because it did demonstrate slightly higher overall velo- 
cities on the beam and also falls nearly into mid-range of 
the cases if D = 15.0 and 20.0 in. are regarded as unreason- 
able models due to their size. These two cases were looked 
at simply to check the effect of the large mass they would 
represent on the hull. Noting Figure 13 that shows the 
maximum velocities of the bracket connection points, the 
rings at frames 18 and 27, it is seen that the variation 
in peak responses is only slight over the range of D = 
4.0-10.0 in. for the tapered charge case. This corresponds 
to substructure to shell mass ratios of . 07644- . 19875 . The 
empty shell velocity at this point is 10.48 ft/sec and its 
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highest maximum velocity occurs when D = 7.0 in. which is 
11.88 ft/sec corresponding to a mass ratio of .13638. The 
conventional charge empty shell velocity at this point was 
32.14 ft/sec and this was the maximum velocity recorded on 
frames 18 and 27 of the conventional cases. 

The velocity-time histories of the beam and bracket for 
the cases D = 2.0-10.0 in. for both charge types are very 
similar to the D = 4.2 in. cases shown in Figures 14 through 
16: At D = 1.0 in. the response was very oscillatory and 

produced the largest maximum velocity for the conventional 
charge at the beam tip. This could be due to the large 
impulse combined with the smallest mass ratio of .01874. 

The other cases simply show overall decreases in velocities 
but retained the basic shape of the veloc i ty-t ime histories. 
Comparing the empty shell to the shell with substructure 
velocity-time histories at the midsection of the shell and 
frame 18 demonstrates what little contribution the 
substructures' responses made in the shell response. In 
Figures 17 through 20 the shape of the veloc i ty-t ime 
history profiles for these points can be seen to be mostly 
unchanged by the presence of the substructure. This was 
true for all cases, conventional and tapered, and increasing 
the weight of the substructure served mostly to reduce the 
velocity at the connection point. 

Comparing the frequencies of the shell model from Table 
II with the first five natural frequencies of the substruc- 
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ture models there are very few matched frequencies at and 
below 600 Hz which are the range expected to excite the 
first two modes of the substructure. The frequencies for 
the cases listed in Table III show how quickly they increase 
to magnitudes where their contributions to response is 
minimal. Any modes of the substructure being excited by the 
base input at the connection points of the brackets and 
shell did not result in a dynamically amplified response. 
Each increment of substructure size produces a set of fixed 
base frequencies of larger magnitude than the previous set. 
If the natural frequency is equal to: 

« = V K/M (7) 

. n e 

where K = stiffness and = effective mass then the 

e 

stiffness must be increasing more rapidly than the mass of 
each structure to produce the higher frequencies. 

B. OVERALL SHELL RESPONSES 

The shell response overall seems to exhibit some of the 
expected modes for the circumferential waves used in the 
calculations with the exception of the purely torsional 
modes that were deleted early on. First looking at the 
velocity-time histories in Figures 21 through 24 the extreme 
difference in velocities produced by the two types of 
charges is apparent. As expected, the conventional charge 
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produces larger velocities in each case. These velocities 
are in the athwart ship direction, positive sense to port, 
at call out points on the midships station and the three 
frames where the discrete rings are located. The velocity- 
time histories due to the conventional charge, illustrate how 
the velocities propagate from the midships station to the 
end plates. Looking at the peak velocities the shell can be 
seen to be bending about the midships station while the 
entire shell translates in the • di rect ion of the shock wave. 
These modes can be seen in the velocity profiles for 
conventional and tapered charge cases where the midships 
peak velocity is less than those for the two rings located 
at frames 9 and 18. Looking at Figures 22 and 25 shows that 
the velocities at the midships station for the conventional 
case on the port and starboard sides are of opposite signs 
at about time equal to 1. 0-2.0 msec. This illustrates a 
higher order circumferential wave of N = 2. Figures 25 and 
26 show the velocity profiles at the top of the shell in the 
athwart ship direction. In each case run the top and bottom 
velocity profiles were exactly matched indicating the 
symmetric response to the shock loading in the translational 
and whipping modes. Figures 27 and 28 illustrate the slight 
lag in the velocities on the port side. As a rule the 
conventional starboard side velocities at ring frames were 
always higher than the port. This is due to the short 
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duration of shock front of the conventional charge. 
Alternately, the tapered charge caused slightly higher peak 
velocities on the port side but also lagged the peak 
responses on the stbd side. This shows how the tapered 
charge allows the port side to exhibit more response to the 
exciting forces. 

The end plate velocities are shown in Figures 29 and 30. 
The end plate in Figure 29 can be seen to bulge in and out 
at the center for the conventional case and in Figure 30 the 
tapered case shows the same response. The ring at frame 
zero tends to hold the outer circumference of the end plate 
in place as shown by the lower velocities experienced at the 
outer edge. These velocity-time histories are indicative of 
breathing modes of the shell. 

C. FORCES IN THE SUPPORTS 

The axial forces were estimated for the D = 4.2 in. 
model subjected to both the conventional and tapered 
charges. In both cases the forces fluctuate between 
compressive and tensile as indicated by the changes in sign 
in the results in Appendix C. The maximum force estimated 
occurred in the conventional case. It was a compressive 
force of about 103 kips or a stress of about 16 psi . The 
maximum tensile force experienced by the support bracket for 
this case was estimated at approximately 75 kips. The 
tapered charge case produced much lower force estimates. 
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The maximum was tensile and occurred at about 1.7 msec 



whereas the conventional maximum force occurred at about 
.3 msec. The maximum estimated tensile and compressive 
forces for the tapered charge case are 25 kips and 23 kips, 
respectively . 
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TABLE IV 



SUMMARY OF RESULTS — TAPERED CHARGE 



D ( i n ) Fl(Hz) F2(Hz) 



MAXIMUM VELOCITIES (FPS) 



Ml /M 2 









TIP( 04 ) 


MID( 21 ) 


J CN (12) 


2/21 (180°) 


1 


25.46 


81 .63 


27.92 


19.76 


10.96 


10.32 


.01874 


2 


50.92 


162.2 


24.58 


21 .79 


11 .39 


10.19 


.0377 


3 


76.38 


240.6 


22.10 


19.14 


11 .02 


10.11 


.05696 


4 


101.8 


315.9 


21.18 


18.85 


11 .94 


11 .37 


.07644 


4.2 


106.9 


330.5 


20.71 


18.52 


1 1 .97 


1 1 .53 


.0765 


5 


127.3 


387.1 


18.98 


17.21 


11 .98 


11.76 


.09602 


6 


152.8 


453.4 


17.49 


15.87 


11 .89 


1 1 .72 


.11615 


7 


178.2 


51 4.0 


15.69 


14.33 


12.12 


11 .88 


.13638 


8 


203.7 


568.7 


1 3.72 


12.70 


11 .83 


11 .69 


.15685 


9 • 


228.7 


616 .3' 


12.05 


11 .76' 


1 1 .43 


11 .24 


.1775 


1 0 


254.6 


659.6 


12.04 


11 .75 


11.11 


10 .9Q 


.19875 


1 5 


381 .9 


798.5 


11.16 


11.17 


10.74 


10.81 


.3071 


20 


509.2 


863.3 


10.79 


10.75 


10.68 


10.37 


.4219 



(Shell mass M2 = 3420.0 lbm. 

Empty shell V(2/21) = 10.48 FPS) 



55 



TABLE V 



SUMMARY OF RESULTS — CONVENTIONAL CHARGES 



D ( i n ) Fl (Hz) 


F2 ( Hz) 


MAXIMUM VELOCITIES (FPS) M1/M2 

TIP ( 0 4 ) M I D ( 2 1 ) JCN (12) 2/21(180°) 


1 


25.46 


81 .63 


77.21 


43.93 


33.97 


31 .75 


.01874 


2 


50.92 


162.2 


55.66 


41 .32 


31 .97 


29.81 


.0377 


3 


76.38 


240.6 


51 .59 


39.90 


33.44 


27.85 


.05696 


4 


101 .8 


315.9 


58 .51 


48.48 


30.77 


28 .02 


.07644 


4.2 


106.9 


330.5 


59.16 


48.91 


29.87 


27.83 


.0765 


5 


127.3 


387.1 


58.04 


48.76 


26.56 


26.43 


.09602 


6 


152.8 


453.4 


54.42 


47.06 


24.18 


23.45 


.11615 


7 


178.2 


514.0 


49.22 


43.05 


24.48 


21 .78 


.13638 


8 


20 3.7 


568.7 


43.33 


38.35 


25.06 


23.02 


.15685 


9 


228.7- 


616.3 


37.78 


33.81 


34.55 


23.33 


.1775 


10 


254.6 


659.6 


32.35 


29 .65 


23.48 


22.90 


.19875 


1 5 


381 .9 


798.5 


20.91 ' 


20.34 


19.65 


19 . 69 


.3071 


20 


509.2 


863.3 


17.44 


17.47 


17.12 


16.95 


.4219 



(Shell mass M2 = 3420.0 lbm, 

Empty shell V(2/21) = 34.14 FPS) 
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MAX VELOCITIES VS BEAM DEPTH 
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Figure 10. Maximum Velocities at Tip of Beam 
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Figure 11. Maximum Velocities at Midsection of Beam 




MAX VELOCITIES VS BEAM DEPTH 
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Figure 12. Maximum Velocities at Beam/Bracket Junction 
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Figure 13. Maximum Velocities at Shell/Bracket Connection at Frame IB 
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Figure 14 . 



Velocity-Time Histories at Beam Tip for D = 4.2 
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Figure 15. Velocity-Time Histories at Beam Midsection for 
D = 4.2 In. 
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VELOCITY PROFILE 



JCNU2) V-flTH 




Figure 16. Velocity-Time Histories at Beam/Bracket Junction 
for D = 4.2 In. 
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VELOCITY PROFILE 



NODE 2/21 V-flTH PORT 




Figure 17. Velocity-Time Histories at Frame 18 for 
D = 4.2 In. 
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VELOCITY PROFILE 



NODE 2/21 V-flTH PORT 




Figure 18. Velocity-Time Histories at Frame 18 for Emoty 
Shell 
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VELOCITY PROFILE 



NOOE 2/21 PORT V-flTH 




Figure 19. 



Velocity-Time Histories at Frame 18 for 
Conventional Charge 
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VELOCITY PROFILE 



NODE 2/21 PORT V-RTH 




Figure 20. Velocity-Time Histories at Frame 18 for 
Tapered Charge 
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VELOCITY PROFILE 



NODE 2/27 V-flTH PORT 




Figure 21. Velocity-Time Histories at Midships with 
Substructure of D = 4.2 In. 
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VELOCITY PROFILE 



NODE 2/27 V-flTH PORT 




Figure 22. Velocity-Time Histories at Midships for Emoty 
Shell 
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VELOCITY PROFILE 



NOOE 2/12 V-RTH PORT 




Figure 23. Velocity-Time Histories at Frame 9 for Emoty 
Shell 
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VELOCITY PROFILE 



NOOE 2/03 V-flTH PORT 




Figure 24. Velocity-Time Histories at Frame 0 for Empty 
Shell 
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V'ELOC I TY PROFILE 



NOOE 2/27 V-flTH CONV 




Figure 25. 



Velocity-Time Histories at Midships for Empty 
Shell from Conventional Charge 
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VELOCITY PROFILE 



NOOE 2/27 V-RTH TAP 




Figure 26. Velocity-Time Histories at Midships for Empty 
Shell from Tapered Charge 
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VELOCITY PROFILE 



NODE 2/21 V-RTH CONV 




Velocity-Time Histories at Frame 18 of Empty 
Shell from Conventional Charge 
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Figure 27. 
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Figure 28. Velocity-Time Histories at Frame 18 of Empty 
Shell from Tapered Charge 
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VELOCITY PROFILE 



END PLATE CONV V-FWD 




Figure 29. Velocity-Time Histories at End Plate for 
Conventional Charge 
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VELOCITY PROFILE 



END PLATE TAP V-FWO 




Figure 30. Velocity-Tine Histories at End Plate for Tapered 
Charge 
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VI. CONCLUSIONS 



The results from the cases studied tend to show that a 
model can be designed to minimize dynamic amplification. 

In no case did a truly resonant situation occur. As the 
substructure model increased in size and mass so did the 
stiffness and the effect was as if a lumped mass was 
attached to the shell. Even though substructure sizes 
having frequencies that matched those of the hull were 
tested, the corresponding mode shapes being excited did not 
result in dynamic amplification of the beam responses. The 
beam model having a secondary frequency of 330.5 Hz became 
slightly amplified but only for the conventional explosive 
analys i s . 

Though the models used in this study are based on real 
configurations it is important to keep in mind this was a 
numerical study on a simplified model of an already simple 
structure. Even though the substructure design used did not 
become highly excited when combined with the shell in this 
study, the beam's behavior might change if combined with a 
shell that has a frequency range and mode shapes that are 
different from those of the shell used in this study. The 
shell model was a relatively small structure that resulted 
in a majority of frequencies too high to excite the 
substructure significantly. 
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ELSHOK has the capability of representing much more 
complex models than those used here. It has the capability 
of simulating a full sized submarine hull with internally 
attached equipment. The ELSHOK code has been shown to give 
fairly accurate results, especially in the early portion of 
shock problems, through comparisons of its predictions to 
physical test results. However, like any numerical device 
it needs such comparisons to verify its results. On the 
other hand, conducting tests using ELSHOK has the distinct 
advantage of eliminating risk of damage to valuable assets. 
Once a model has been verified to satisfactorily predict the 
prototype responses, a variety of internal equipment models 
could be applied to perform a variety of studies at 
fractions of the cost of building real models. Again the 
drawbacks are the pitfalls of any numerical endeavor, 
namely, the output is only as good as the input. Simplified 
models may produce simplified results and poor assumptions 
must be guarded against. 

It is important to realize that the results of this 
study were arrived at through the implementation of a 
potentially powerful tool for future research at the Naval 
Postgraduate School, namely ELSHOK. Further studies in the 
area of UNDEX using the code should be carried out to prove 
its validity or develop it into more useful applications of 
design and research. Suggestions on improvements to develop 
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the line of research in this thesi 
sensitivity study into a realistic 
construct a representative full-si 
for B0S0R4 and model all internal 
appropriate detail while applying 
degrees of freedom. 



s from a numerical 
problem are to obtain or 
ze submarine model or SSTV 
equipments to be tested in 
fewer restrictions to the 
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APPENDIX A 



INPUT FILES FOR ELSHOK 

1. BOSOR4 INPUT DATA FOR ALL CASES 

The following data file was used to generate the 
cylindrical shell model. Four separate runs of B0S0R4 were 
made corresponding to N = 0,1, 2,3 by setting NOB, NMINB and 
NMAXB equal to each case N. All other inputs remain the 
same for each case. 
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RING -STIFFENED 
2 



H 

H 



1 

0 

2 



7 
3 
1 
0 
0 

16. 81250 
0 
0 
3 

2. 500000 
N 

0 

0 

0 

0 

2 

0. 1080000E+08 

0. 3200000 

0. 2535000E-03 
0 
0 
1 
N 
N 
H 
H 

25 

3 

1 

16. 81250 
0 

16. 81250 
129. 0625 
0 
3 

0. 1875000 
N 

3 

2 

3. 500000 
53. 12500 
103. 7500 
1 
1 
1 

0. 3000000E+08 
7. 125000 
7. 535500 
2. 375000 
0 

-1. 781250 
0 

. 7501E+08 

0. 7330000E-03 

0. 3000000E+08 
8. 250000 
11. 69800 
2. 750000 
0 

-2. 062500 
0 

0. 9249000E+08 

0. 7330000E-03 

0. 3000000E+08 
8. 250000 



^INDER WITH FLAT ENDPLATES-- SYMMETRY USED--N=0 
INDIC = ANALYSIS TYPE INDICATOR 

NPRT = OUTPUT OPTIONS (1=MINIMUM, 2=MEDIUM, 3=MAXIMUMJ 
ISTRES= OUTPUT CONTROL ( 0=RESULTANTS , l=SldMA, 2=EPSILON) 
NSEG = NUMBER OF SHELL SEGMENTS (LESS THAN 25) 



SEGMENT NUMBER 



NMESH = NUMBER OF NODE PQINTS (5 = MIN. , ,, 

1 . OR 2 OR 3) F6R NODAL POINT SPACING 



NTYPEH= CONTROL INTEGER 



111 
98 = MAX. )( 1) 



NSHAPE= INDICATOR (1,2 OR 4) FOR GEOMETRY OF MERIDIAN 
R1 = RADIUS AT BEGINNING OF SEGMENT f SEE P. 66) 

Z1 = AXIAL COORDINATE AT BEGINNING OF SEGMENT 

R2 = RADIUS AT END OF SEGMENT 

Z2 = AXIAL COORDINATE AT END OF SEGMENT 

IMP = INDICATOR FOR IMPERFECTION (0=NONE, l=SOME) 

NTYPEZ = CONTROL ( 1 OR 3 ) FOR REFERENCE SURFACE LOCATION 
ZVAL = DISTANCE FROM LEFTMOST SURF. TO REFERENCE SURF. 

DO YOU WANT TO PRINT OUT R(SJ. R T (S), ETC. FOR THIS SEGMENT? 
NRINGS= NUMBER ?MAX=20) OF DISCRETE AlNGS IN THIS SEGMENT 
K=ELASTIC FOUNDATION MODULUS (E.G. LB/IN**3)IN THIS SEG. 
LINTYP= INDICATOR (0, 1, 2 OR 3) FOR TYPE OF LINE LOADS 
NLTYPE=CONTROL (0,1, 2, 3) FOR TYPE OF SURFACE LOADING 
NWALL= INDEX (1, 2, 4, 5, 6, 7, 8) FOR WALL CONSTRUCTION 

e = young m6duAus foA sAin 

U = POISSON'S RATIO FOR SKIN 

SM =MASS DENSITY- OF SKIN (.E.G. ALUM. =.00025 LB-SEC**2/IN**4) 
ALPHA = COEFFICIENT OF THERMAL EXPANSION 

ANRS = CONTROL fO OR 1). FOR ADDITION OF SMEARED STIFFENERS 
SUR = CONTROL FOR THICKNESS INPUT (0 OR 1 OR -1) 

DO YOU WANT TO PRINT OUT THE C(I,J) AT MERIDIONAL STATIONS? 
DO YOU WANT TO PRINT OUT DISTRIBUTED LOADS ALONG MERIDIAN? 



SEGMENT 

NMESH = 

NTYPEH= 

NSHAPE= 

R1 

Z1 

R2 

Z2 

IMP = 
NTYPEZ = 
ZVAL = 
DO YOU 
NRINGS= 
NTYPE = 
Z( I) = 
Z I = 
Zflj = 
NTYPER= 
NTYPER= 
NTYPER= 
E = 

A = 

IY 

IX = 

IXY 

El 

E2 = 

GJ 

RM 

E = 

A = 

IY 
IX 

IXY = 

El 

E2 = 

GJ 

RM = 

E = 

A = 



NUMBER 2 2 

NUMBER OF NODE P 



2 

)INTS 



2 

(5 = MIN. ; 98 = MAX. )( 2) 

CONTROL INTEGER ( 1 OR 2 OR 3 ) f6r NODAL POINT SPACING 



INDICATOR (1,2 OR 4) FOR GEOMETRY OF MERIDIAN 
= RADIUS AT BEGINNING OF SEGMENT (SEE P. 66) 

= AXIAL COORDINATE AT BEGINNING OF SEGMENT 
= RADIUS AT END OF SEGMENT 
= AXIAL COORDINATE AT END OF SEGMENT 
INDICATOR FOR IMPERFECTION (0=NONE, l=SOME) 

CONTROL ( 1 OR 3 ) FOR REFERENCE SURFACE LOCATION 
DISTANCE FROM LEFTMOST SURF. TO REFERENCE SURF. 

WANT TO PRINT OUT R( S ) , R T (S), ETC. FOR THIS SEGMENT? 
NUMBER (MAX=202 OF DISCRETE &INGS IN THIS SEGMENT 
CONTROL FOR IDENTIFICATION OF RING LOCATION (2=Z, 
AXIAL COORDINATE OF ITH RING, Z( 1) 

AXIAL COORDINATE OF ITH RING, Z ( 2) 

AXIAL COORDINATE OF ITH RING, Z( 3] 

TYPE (0 OR 1 OR 2 OR 4 OR 5) OF DISCRETE RING NO. 
TYPE (0 OR 1 OR 2 OR 4 OR 5 1 

TYPE (0 OR 1 OR 2 OR 4 OR ' 
young's MODULUS OF RING? 1 
CROSS SECTION AREA OF RING 

MOMENT OF INERTIA ABOUT Y- 

MOMENT OF INERTIA ABOUT X-AXIS( 

PRODUCT OF INERTIAL l] 

RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 

AXIAL COMPONENT OF RING ECCENTRICITY( 1) 

TORSIONAL RIGIDITY? 1) 

RING MATERIAL DENSITY (E. G. ALUMINUM= . 0002535 ) ( 1) 
YOUNG'S MODULUS OF RING( 2) 

CROSS SECTION AREA OF RING? 2) 

MOMENT OF INERTIA ABOUT Y-AXIS (SEE FIG. 

MOMENT OF INERTIA ABOUT X-AXIS( 2) 



3 = R) 



OF DISCRETE RING NO. , 
OF DISCRETE RING NO. ( 

CI^ ( SEE FIG. ON P. 70 )( 
1) 



1) 



70)( 1) 



ON P. 70 ) ( 2) 



PRODUCT OF INERTIA( 2) 

RADIAL COMPONENT OF RING ECCENTRICITY (SEE P. 
v 2) 



70 )( 2) 



AXIAL COMPONENT OF RING ECCENTRICITY( 

TORSIONAL RIGIDITY? 2), w „ 

RING MATERIAL DENSITY ( E. G. ALUMINUM=. 0002535 )( 2) 
YOUNG'S MODULUS OF RING?” 3) 

CROSS SECTION AREA OF RING? 3) 
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? 



;? 



DO YOU WANT TO LIST OUTPUT FOR SEGMENT ( 2) 

DO YOU WANT TO LIST PREBUCKLING RESULTANTS AND RING FORCES? 

DO YOU WANT TO LIST FORCES IN THE DISCRETE RINGS, IF ANY? 
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2. ACESNID INPUT DATA FOR ALL CASES 



The input data for ACESNID must encompass all waves N 
used for the BOSOR4 model. The waves are specified by 
NSTART-NFINISH . 



VIRTUAL MASS FOR RING STIFFENED CYLINDER WITH FLAT ENDS, N=0, 1,2,3 



0 3 1 
12.61 



1 0 /NSTART , NFINISH , NFREQ , NVMASS , NCHECK 
■ . _ 12 2 1 /NFLAT,NCYL,tfRITE,NWBOSG,tfSYMF 
2 3 12 0 0 /NORDER, ftFEND& .NFCE^T, NFCM^T , NOMIT 
9.59684E-5 5. 833E4 0.005 ?RHOFL,^SOUND , £RR 



1111110111 
1.0 /CPSflj 
1119 /SETS JSGBEi 



/OUTPUT FLAGS' 



_ ✓ JSGBEG .JPTBEG , JSGEND .JPTEND LEFT 

2 1 2 27 /SEFS JSGBEG .JPTBEG .JSGEND , JPTEND CENTRAL 
1119 /SOURCES JSGBEG, JPTBEG, JSGEND, JPTEND LEFT 
2 1 2 27 /SOURCES JSGBE& , JPTBE& , JSGENf) , JPTEND CYLINDER 
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3. PIFLASH INPUT DATA FOR ALL CASES 

BOSSY100-103 are the shell mode files produced from the 
BOSOR4 runs. Note the deleted shell modes that were purely 
torsional for N = 0. 



4 0 0 1 
2 2 2 2 
386.4 
B0SSY100 
BOSSY101 
BOSSY102 
BOSSY103 
22 30 * 

3 i l 

26*27 
12 3 
26 27 
12 3 
26 27 28 



/NUMBER NTORSN NPTM NSYM NSYMP 
/(NWETSEG(K) ,NUSESEG(K) ,K=1,2 

/Gravity 

/ SMF ,N=0 
/SMF N=1 
/SMF N=2 
7 SMF !n=3 

30 30 / (NJUSEfJ) ,J= l.NITEMS) 

4 / (korGs(k),k=1’nko^gs) 

7 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 

' ‘ - ----- u 12 13 14 - ------ -- -- 



4 5 6 7 8 9 10 
28 29 30 
4 5 6 7 8 9 10 
28 29 30 
4 5 6 7 8 9 



30 



11 12 13 14 15 



- . - .10 11 12 13 14 15 
29 30 / (JUSE(J,K), J=NJUSE 



16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


16 


17 


18 


19 


20 


21 


22 


23 


24 


25 


.1.6 


17 


18 


19 


20 


21 


22 


23 


24 


25 
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4. SAP IV INPUT DATA FOR D = 4.2 IN. CASE 

The following data file generated the internal structure 
model. Note that only half the doubly overhung beam is 
modeled as indicated by the beam lengths at the node, i.e., 
node 21 -- 44.81250. 



SYMMETRICAL SUBSTRUCTURE- -FOR BEAM OF D=4. 2 IN 
21.1.0. 10 100 0.0.0 /NUMP,NELTYP. . . NRIGID 



0 

0 

0 

0 

0 

0 

0 

2 

1 

2 
1 



1 1 1 0 1 1 1 20.0 0.0 16.71875 
3 1 I 0 1 1 1 20.0 0.0 5. 5729167 
4110101 0. 0 0. 0 0.0 

1 20.0 0.0 0 . 0 

1 22. 756944 0. 0 0. 0 
1 42. 055556 0. 0 0.0 
1 44.81250 0.0 0.0 

BEAMS- 20 ELEMTENTS TOTAL 
7.29E-04 0.0/BEAM MOD, POIS , DENS 
7. 29E- 04 0.0 /BRACKET " 



12 

13 

20 

21 

20 



110 
110 
110 
110 
,.202 
29. 0E+06 
29. 0E+06 



1 0 
1 0 
1 0 
1 1 

U 

3 



0.0 /GEOMETRY AND DOF 
1 0. 0 
0.0 
0. 0 
0. 0 
0. 0 
0. 0 



8.4 0.0 0.0 1.0 2.8 12.35/BEAM AREA AND MOMENTS OF I 
6.3 _0._0_0.0_1. 0 20.84 .525/BRACKET ir w 



2 ....... 

0. 0 0. 0 0. 0 0. 0 

0. 0 0. 0 0. 0 0. 0 

0. 0 0.0 0.0 0.0 / MUULT G-X,Y,Z 
' ' 2 21 2 2 0 0 0 0 0 0 0 
3 21 220000000 

12 21 2 2 0 0 0 0 000000 000001 0/SET NODE RELEASE AT 12 
51110000000 
20 20 21 1110000000 
0 0 0.0 0.0 0.0 0.0 0.0 0.0 /CONC MASSES 
0. 0 0. 0 0. 0 0. 0 /ELEMENT LOAD MULTIPLIERS 

0 0 0 0. 0 0. 0 0 0. 0 /IFPRjIFSS,. . . ,NFO, SHIFT 

1 1 /NFACE.NSYM 

1 1 1 1 1 1 1 /NODE , ( FIXITYf J ) , J= 1, 6 
END DATA FOR D (DUMMY TITLES) 

0000000 000 



2 2 

u 
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5. 



PICRUST INPUT DATA FOR D 



4.2 IN. CASE 



The circumferential location of the substructure 
attachment point is specified by ANDEG which is the angle 
ccw from the shell global Z-axis. The meridional location 
is specified by the node 21 at frame 18. The angle DEGROT 
specifies the orientation of the substructure local z-axis 
to the shell Z-axis. This is indicated by in Figure 5. 



1111011100 

0111110101 /OUTPUT FLAGS 1-20 
500 /NJUSE , NHWS0B , NHWBAR 
1 2 21 180.0 /NIP SUB, LB0SEG ,LB0SPT , ANDEG 
180.0 / DEGROT 

12345/ ( JUSE( J) ,J=1, NJUSE) 
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6. USLOB INPUT DATA 

In Phase III of ELSHOK the USLOB code perforins the 
integration of the incident pressures using a Runge-Kutta 
scheme. NCHARGE specifies which type of UNDEX is to be 
represented. For the tapered charge nine discrete pressure- 
time points are used as input. The numbers following the 
charge weight of 352 lb. are the constants used in equations 
(1) and (2) that empirically describe the explosive HBX-1 
(order-Kl, K2, K3, K4) . Note the 0.0 for surface cutoff 
effects. For the empty shells NSUBS must be set to 0. The 
integration is performed to compute velocities at all the 
nodes that make up the shell and substructure but only -a 
couple are requested as printed output in NPTSH and NPTSUB 
as indicators that the program ran successfully. Nodes 1/01 
and 2/21 are checked on the shell and node 1 and 21 are 
checked on the substructure. The value NQUAN = 3 requests 
the z direction velocities as output. 



TA 331 E 2 2"4 R i E l 1 /NTIME .NSKIP .NCHARGE .NQUAD.NFINE.KOUPLE, NSUBS 
3:0E-05 129.0625 856. 8l25 0.6 /DELT,XLUAD,RLOAD,SURCUT 

9 /NSHAPE „ 00053 350.0 .00243 305.0 .00294 160.0 

00540 48.0 i 00565 0.0/DISCRETE POINT INPUTS 



0 . b 240. 0 . 0002 320. 0 
.00340 115.0 .0044 70.0 
0.0 /DECAY 

1 
2 
1 
2 



1010 /OUTPUT FLAGS 1-20 



0 10 0 1 
0 10 10 
/NPTSHL 
1 2 21 /LBOSEG , LBOSPT 
/NPTSUB 



1 3 21 3 /NSTAT, NQUAN 
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CONVENTIONAL CHARGE 

331 2 1 9 5 1 1/NTIME.NSKIP.NCHRG,. . . ,NSUBS 
3.0E-05 129.0625 856. §125 0.0/ DELT , XLOAD , RLOAD , SURCUT 
352 3. 835E+05 1.144 3.031E-05 -0.247 
1010101011 
111000101 1/OUTPUT FLAGS 
2/NPTSHL 

112 2 1 / LBO S EG , LBSPOT 
2/NPTSUB/ NSTAT , 1&QUAN 
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7. PUSLOB INPUT DATA 

This data file generates the velocity profiles that are 
plotted on Tektronix devices. The time step DELT and the 
interval of integration specified by NTIME and NSKIP must be 
the same as for USLOB. This input requests velocity 
profiles be expressed in terms of the shell global 
coordinate system. Shell output is requested by segment 
node and angle. Each substructure point is specified by 
node, number of degrees of freedom at the node for which 
velocity profiles are to be produced and which degree of 
freedom it is produced for.' The tip is node 4, one degree 
of freedom in the local z direction. 



331 211 l/NTIM, NSKIP, NSUBS , NTEK , NCARD 

5iM? E Sl»/s E c. 

16 / NPTSHL-W 
V-ATH 1/01 
11-11 0.0 
V-ATH 1/04 

H -n 0.0 180.0 
V-ATH 1/09 „ 

1 9 1 2 0. 0 180. 0 

V-ATH 2/03 „ 

2 3 1 2 0. 0 180. 0 
V-ATH 2/07 

2 7 1 2 0.0 180. 0 
V-ATH 2/12 „ „ 

2 12 1 2 0. 0 180. 0 
V-ATH 2/16 „ „ 

2 16 1 2 0.0 180. 0 
V-ATH 2/21 
2 21 1 2 0. 0 180. 0 
V-ATH 2/27 „ 

2 27 1 2 0. 0 180. 0 
V-ATH 1/09 

1 9 1 2 90.0 270.0 
V-ATH 2/03 ^ 

2 3 1 1 90.0 
V-ATH 2/07 

2 7 1 1 90.0 
V-ATH 2/12 
2 12 1 1 90.0 
V-ATH 2/16 
2 16 1 1 90.0 
V-ATH 2/21 
2 21 1 1 90. 0 

2"27 H l 2 { 2 9o! f O U /LBOSEG,LBOSPT^NSPHC,NANG,ANGDEG 
16/NPTSHL-U 
V-FWD 1/01 
1 1 1 1 0.0 
V-FWD 1/04 
1 4 1 2 0. 0 180. 0 
V-FWD 1/09 
1 9 1 2 0. 0 180. 0 
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V-FWD 2/03 
2 3 -1 2 0.0 180.0 
V-FWD 2/07 
2 7 -1 2 0. 0 180. 0 
V-FWD 2/12 
2 12 -1 2 0.0 180.0 
V-FWD 2/16 
2 16 -1 2 0. 0 180. 0 
V-FWD 2/21 
2 21 -1 2 0.0 180.0 
V-FWD 2/27 
2 27 -1 2 0.0 180.0 
V-FWD 1/09 

1 9 1 2 90. 0 270.0 
V-FWD 2/03 

2 3 -1 1 90.0 
V-FWD 2/07 

2 7 -1 1 90.0 
V-FWD 2/12 
2 12 -1 1 90.0 
V-FWD 2/16 
2 16 -1 1 90.0 
V-FWD 2/21 
2 21 -1 1 90.0 
V-FWD 2/27 
2 27 -1 1 90. 0 
16/NPTSHL-V 



V-DWN 1/01 
1 1 -1 1 0.0 
V-DWN 1/04 
14-12 0.0 180.0 
V-DWN 1/09 

1 9 -12 0. 0 180. 0 
V-DWN 2/03 

2 3 -1 2 0.0 180.0 
V-DWN 2/07 

2 7 -1 2 0.0 180.0 
V-DWN 2/12 
2 12 -1 2 0.0 180.0 
V-DWN 2/16 
2 16 -1 2 0.0 180.0 
V-DWN 2/21 
2 21 -1 2 0.0 180.0 
V-DWN 2/27 
2 27 -1 2 0.0 180.0 
V-DWN 1/09 

1 9 -1 2 90. 0 270.0 
V-DWN 2/03 

2 3 -1 1 90.0 
V-DWN 2/07 

2 7 -1 1 90.0 

V-DWN 2/12 

2 12 -1 1 90.0 

V-DWN 2/16 

2 16 -1 1 90.0 

V-DWN 2/21 

2 21 -1 1 90.0 

V-DWN 2/27 

2 27 -1 1 90. 0 

3/NPTSUB 

TIP ( 04J/NGAGE 

jcn 1 / J^STAT.NDOFPT.NQUAN 

12 13 



MID (^21) 



21 



93 



nnnnn 



APPENDIX B 



PROGRAMS USED TO PREPARE AND PLOT VELOCITY PROFILES 



PROGRAM PLOTCONV 

THE PURPOSE OF THIS PROGRAM IS TO READ THE PUNCH-CARD 
FILE FROM PUSLOB AND CONVERT IT TO A FORMAT WHICH CAN 
BE ACCEPTED BY EASYPLOT OR DISSPLA ON IBM 3033. 

DIMENSION VEL( 166 , 32 ) , VELO(166,2), LABEL( 32) , TAG(32), LTITLE( 32 ) 

CHARACTER'* 10 LABEL 

CHARACTER*40 LTITLE 

COMMON / CSOLVE / DTRECS 

NCURV=32 





LL=L*20 
DO 30 M 
VELO(M, 



END 
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PROGRAM FISH1 IS THE FORTRAN CODE REQUIRED TO PLOT 

VELOCITY PROFILES USING DISSPLA 

REAL*4 VEL( 166 ) , TIM( 166 ) , TIM2(166 ) , VEL2( 166 ) 

DO 10 1=1,166 
READ(20,*)tIM(I),VEL(I) 

10 CONTINUE 

DO 20 J= 1 , 166 

READ(_21 , * )TIM2( J ) ,VEL2( J) 

20 CONTINUE 
CALL COMPRS 
CALL PAGEf 8. 0,8. 0) 

CALL AREAZD(6. ,5. $) 



CALL XNAMEt 'TIME (MSEC)?’ 100] 

CALL YNAME? VELOCITY ( FT/$Ecl? ’ , 100 ) 

CALL HEADIN( 'VELOCITY PROFILE ?', 100 , 2. , 



c ,f. 0. 



( S EMPTY SHELLS ’ , IPAK, 1) 
('LOADED SHELL? ' , IP a£ , 2 

1] 

' $’, 100 ) 



CALL FRAME 
CALL INTAXS 
CALL YAXANG( 0 
CALL YTICKS " 

CALL XTICKS s _ a 
CALL GRAF( 0. , 2 . 0,10. .-10. ,10. ,35. ) 
CALL MESSAG( 'NODE 2/21 PO&T V-ATH 
CALL LINESP ( 1. 75) 

CALL LINES 

CALL LINES 

XW=XLEGND ( tPAK , 2 
YW=YLEGND ( IPAK.2 
CALL MYLGEN ( ’KliY 
CALL LEGLIN 
CALL POLY 3 
CALL DASH 
CALL CURVEf TIM, VE: 

CALL RESETl t DASH T . 

CALL CURVE(TIM2,V£L2, 166,0) 

CALL DOT 

CALL GRID( 1 , 1 ) , 4 
CALL RESExf 'DOT ) 

CALL LEGEND (IPAK , 2 , 4. 0 , 4. 5 ) 

CALL ENDPLiO ) 

CALL DONEPL 
STOP 
END 



1 ) 



C$’ ,23,1. ,5. 7) 



) 



.,166,0) 
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APPENDIX C 



PROGRAM USED TO ESTIMATE AXIAL FORCES IN THE 

SUPPORT BRACKETS 



PROGRAM FORCES 

ESTIMATES FORCES DEVELOPED IN THE SUPPORT BRACKETS 
USING SIMPS 1/3 INTGRTR TO EST DEFLECTIONS 
REAL-4 TT( 166 ) , TB(- 166 ) , VT( 166 , 2} , VBU66 , 2 ) , xf 166 , 2 ) ,XDEF( 166 
REAL*4 XL.EMOD,H,AREA,VTO,VTl,VT2,VBO,VBi;VB2,DT DB,DEF 
REAL*4 XFdRC,TiM£,D 



) 



XL= 10. 4 
EMOD=. 3E+08 
D=4. 2 

AREA=D*1. 5 
H=6. 0E-05 
DO 10 1=1,166 

READ(81,*) ( VT( I ,M) ,M= 1 , 2 ) 
READ(82,*) (VB(I,N) N=l,2) 
10 CONTINUE 

DO 20 J= 1 , 165 , 2 
II=J 
JJ=II+1 

kk=ii+; 

' '.II. 2 



VT0=VT 
VT1=VT 
VT2=VT 
VB0=VB 
VB1=VB 
VB2=VB 
DT=( VT 



DB 

DEF 



if 



JJ 2 
KK , 2 
II 2 
JJ 2 
KK 2 

+ 4.’~tfTl+VT2)*H*4. 0 



BO+4. *VB1+VB2 )*H*4. 0 



20 



25 



30 



, DT-DB) 

XFORt=( DEF7XL)*EMOD*AREA 
XLJ i 2)=XFOkC 
TIm6=VT( JJ, 1) 
xfj,l)=TIMt 
XDEkfj)=DEF 
CONTINUE 

DO 25 LL= 1,163,2 

WRITEf 6 , * ) txrtL.MM) ,MM=1, 

WRITeU!*) XDEF(f,L) 

CONTINUE 

3.*) ’ MSEC 



2 ) 



WRITE( 83 
DO 30 L= 
WRITE (8 3 
CONTINUE 
STOP 
END 



,NN ) ,NN= 1,2) 



LBF 



TAPERED CHARGE 



CONVENTIONAL CHARGE 



MSEC 

5. 9999999E-02 
0. 1800000 
0. 3000000 
0. 4200000 
0, 5400000 
0, 6600000 
0, 7800000 
0. 9000000 
1. 020000 
1. 140000 
1. 260000 
1. 380000 
1. 500000 
1. 620000 
1. 740000 
1. 860000 
1. 980000 
2. 100000 
2. 220000 
2. 340000 



LBF 

245. 0312 
-4871. 925 
-19920. 20 
-22457. 08 
-16749. 71 
-5310. 641 
5320. 271 
4192. 220 
7970. 726 
13106. 43 
62. 28399 
-14969. 11 
-13569. 19 
3981. 250 
25546. 64 
24680. 42 
8301. 491 
5590. 140 
9395. 421 
9359. 116 



MSEC 

5. 9999999E-02 
0. 1800000 
0. 3000000 
0. 4200000 
0. 5400000 
0. 6600000 
0. 7800000 
0. 9000000 
1. 020000 
1. 140000 
1. 260000 
1. 380000 
1. 500000 
1. 620000 
1. 740000 
1. 860000 

1. 980000 
2. 100000 
2. 220000 

2. 340000 



LBF 

1461. 115 
-30141. 24 
-102820. 2 
-81457. 16 
-31330. 44 
35335. 74 

74938. 54 

40918. 55 
54392. 00 
65419. 58 

-25939. 18 
-89322. 17 
-47749. 70 
62073. 85 
155272. 9 
95783. 45 
-15545. 27 
-12975. 48 

14247. 56 
13262. 66 



96 



2. 460000 


10169. 69 


2. 460000 


18125. 72 


2 . 580000 


11816. 05 


2. 580000 


22709. 01 


2. 700000 


18961 . 72 


2. 700000 


50505. 12 


2. 820000 


19368. 80 


2. 820000 


21335. 82 


2. 940000 


- 382. 3864 


2. 940000 


- 103091. 1 


3. 060000 


- 13737. 00 


3. 060000 


- 140059.4 


3. 180000 


- 8604. 939 


3. 180000 


- 64635. 81 


3. 300000 


- 2495. 676 


3. 300000 


- 7922. 625 


3.420000 


446. 2275 


3. 420000 


14542. 41 


3. 540000 


- 3972. 923 


3. 540000 


- 10436. 12 


3. 660000 


- 9858. 912 


3. 660000 


- 35627. 13 


3. 780000 


- 5315. 625 


3. 780000 


- 6730. 258 


3. 900000 


- 6726. 142 


3. 900000 


- 28879 . 58 


4. 020000 


- 15112. 29 


4. 020000 


- 71626. 63 


4. 140000 


- 11547. 92 


4. 140000 


- 23481 . 91 


4. 260000 


- 4874. 153 


4 . 260000 


31006. 96 


4. 380000 


- 3058. 313 


4 . 380000 


42378. 32 


4. 500000 


- 4332. 930 


4. 500000 


24937. 62 


4. 620000 


- 11046. 69 


4 . 620000 


- 21828. 76 


4. 740000 


- 8562. 529 


4. 740000 


- 1085. 286 


4. 860000 


2656. 177 


4. 860000 


59473. 24 


4. 980000 


1821. 597 


4. 980000 


39020. 46 


5. 100000 


- 2468. 413 


5. 100000 


13899. 35 


5. 220000 


- 403. 8770 


5. 220000 


27988. 47 


5. 340000 


1906. 643 


5. 340000 


30918. 16 


5. 460000 


6909. 898 


5. 460000 


38729. 76 


5. 580000 


8255. 736 


5. 580000 


17834. 63 


5. 700000 


6907. 107 


5. 700000 


- 8653. 294 


5. 820000 


14599. 11 


5. 820000 


28207. 81 


5. 940000 


16075. 41 


5 . 940000 


29423. 95 


6. 060000 


4580. 446 


6. 060000 


- 29775. 31 


6. 180000 


- 2501.472 


6. 180000 


- 46553.45 


6. 299990 


- 3046. 579 


6. 299990 


- 32686. 68 


6. 419990 


1007. 863 


6. 419990 


- 7771. 345 


6. 539990 


4450. 337 


6. 539990 


5449. 699 


6. 659990 


- 1895. 174 


6. 659990 


- 35164. 12 


6. 779990 


- 4134. 870 


6. 779990 


- 42226. 15 


6. 899990 


592. 2988 


6. 899990 


- 10202. 90 


7. 019990 


- 4900. 276 


7. 019990 


- 33267. 90 


7. 139990 


- 13168. 97 


7. 139990 


- 53038. 88 


7. 259990 


- 14155 . 02 


7. 259990 


- 31495. 45 


7. 379990 


- 10175. 60 


7. 379990 


- 3401. 172 


7. 499990 


- 2037. 099 


7. 499990 


31961. 83 


7. 619990 


- 2537. 893 


7. 619990 


16792. 98 


7. 739990 


- 10365. 33 


7. 739990 


- 23224. 84 


7. 859990 


- 6663. 691 


7. 859990 


10102. 76 


7. 979990 


- 212. 9292 


7. 979990 


42517. 71 


8. 099990 


- 1417. 545 


8. 099990 


26730. 40 


8. 219990 


- 1980. 576 


8. 219990 


21689. 60 


8. 340000 


- 1527. 583 


8. 340000 


. 15827. 63 


8. 460000 


4144. 204 


8. 460000 


31344. 69 


8. 580000 


12648. 72 


8. 580000 


58115. 68 


8. 700000 


9013. 509 


8. 700000 


21346. 57 


8. 820000 


4134. 388 


8. 820000 


- 3069. 043 


8. 940000 


8353. 832 


8. 940000 


24245. 50 


9. 060000 


9588. 924 


9. 060000 


20669. 76 


9. 180000 


8389. 120 


9. 180000 


6704. 038 


9. 300000 


7100. 230 


9. 300000 


- 2779. 436 


9. 420000 


4991. 131 


9. 420000 


- 18157. 43 


9. 540000 


9173 . 441 


9. 540000 


2506. 271 


9. 660000 


9349. 440 


9. 660000 


- 908. 7243 


9. 780010 


- 84. 61229 


9. 780010 


- 47302. 98 



97 



INITIAL DISTRIBUTION LIST 



1. Defense Technical Information Center 
Cameron Station 

Alexandria, Virginia 22304-6145 

2. Library, Code 0142 
Naval Postgraduate School 
Monterey, California 93943-5002 

3. Professor Y. S. Shin, Code 69Sg 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93943-5000 

4. Professor R. E. Newton, Code 69Ne 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93943-5000 

5. Department Chairman, Code 69 
Department of Mechanical Engineering 
Naval Postgraduate School 
Monterey, California 93943-5000 

6. Dr. N. T. Tsai 
Defense Nuclear Agency 
SPSS 

Washington, D.C. 20305-1000 

7. Dr. E. Sevin 
Defense Nuclear Agency 
Washington, D.C. 20305-1000 

8. Dr. M. L. Baron 
Weidlinger Associates 
333 Seventh Avenue 
New York, NY 10001 

9. Dr. Ranlet 
Weidlinger Associates 
333 Seventh Avenue 
New York, NY 10001 



98 



10 . 



1 



Dr. Andrew P. Misovec 
5202 W. Military Hwy 
Chesapeake, Virginia 23321 

11. Dr. B. Whang, Code 1750.2 1 

Hull Group Head, Submarine Protection Div. 

David Taylor Naval Ship Research 
and Development Center 
Bethesda, Maryland 20084 

12. Dr. Huang, Code R14 1 

Naval Surface Weapon Center 

White Oaks 

Silver Spring, Maryland 20910 

13. LT G. F. DeConto 2 

47 Water Street 

Sandwich, Massachusetts 02563 



99 



Thesis 
D234435 
c. 1 



Thesis 
D234435 
c. 1 



21(6 






DeConto 

A parametric study 
of elastic response 
of submar ine- 
installed equipment 
subjected to UNDEX 
side-on loading. 



.'3S •• 

DeConto 

A parametric study 
of elastic response 
of submar ine- 
installed equipment 
subjected to UNDEX 
side-on loading. 



